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Abstract. In 0, we studied p-mean curvature and the associated p-minimal 
surfaces in the Heisenberg group from the viewpoint of PDE and differential 
geometry. In this paper, we look into the problem through the variational 
formulation. We study a generalized p-area and associated (p-) minimizcrs in 
general dimensions. 

We prove the existence and investigate the uniqueness of minimizers. Since 
this is reduced to solving a degenerate elliptic equation, we need to consider the 
effect of the singular set and this requires a careful study. We define the notion 
of weak solution and prove that in a certain Sobolev space, a weak solution is a 
minimizer and vice versa. We also give many interesting examples in dimension 
2. An intriguing point is that, in dimension 2, a C 2 -smooth solution from the 
PDE viewpoint may not be a minimizer. However, this statement is true for 
higher dimensions due to the relative smallness of the size of the singular set. 



1. Introduction and statement of the results 

The p-minimal (or X-minimal, H-minimal in the terminology of some authors, 
e.g., jS], CHI) surfaces have been studied extensively in the framework of geo- 
metric measure theory. Starting from the work j3J, we studied the subject from the 
viewpoint of partial differential equations and that of differential geometry (we use 
the term p-minimal since this is the notion of minimal surfaces in pseudohermitian 
geometry; "p" stands for "pseudohermitian"). 

Let be a bounded domain in R 2n . Let X = (xi, xy, X2, X2>, ■■, x n , x n >) £ f2. 
For a graph (X, u(X)) in the Heisenberg group of dimension 2n+ 1 with prescribed 
p-mean curvature H = H(X), the equation for u : fi C R 2n — > R reads 

(1.1) dw Vu -*: = H 

\Vu-X*\ 

where X* = (xy , —Xi, xy , ~X2, x n i , —x n ) (see (|2.10|) in Section 2 for a geometric 
interpretation). In general, for a vector field G — (31,32, •••,<?2n) on C R 2n , we 
define G* = (32, ~9\, 9i, —33, 32n, — 32n-i)- The equation i|l.l|) is the Euler- 
Lagrange equation (away from the singular set Vm- X* = 0) of the following energy 
functional (called the p-area of the graph defined by u if H — 0, see Section 2): 

(1.2) X(u)= I {|Vit- X*\+ Hu}dx 1 f\dx v A ... A dx n f\ dx n , . 

Jn 
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Since we consider the variation over the whole domain, the singular set will 
cause the main difficulty in the study. In order to explain this, we generalize X(-) 
by considering an arbitrary vector field F — F(X) instead of —X* in the following 
form: 



(1.3) T q {u) = / {|Vlt + F\ q + qHu}dx 1 A dx 2 A ... A dx r 



for 1 < q < co, where ft C R m . Let S(u) denote the singular set of u, consisting of 
the points where Vu + F = 0. Let u £ = u + ecp. It is easy to compute (see Section 
3 for the case q — 1) the first variation of T q : (omitting the Euclidean volume 
element) 



(1.4) 



dT q (u e ) 

=o± 



de 

,[ |V(^|«+ f q\Vu + F\i- 2 (Vu + F)-V<p+ [ qH V 
Js{u) Jn\s(u) Jq 

where c q ~ ±1 for q = 1 and c q = for 1 < q < oo. 

For q = 1, can we ignore the term ± f s , u s |Vy|? A recent paper of Balogh 
answered this question completely. In PP Balogh studied the size of the singular set 
S(u) (called the characteristic set in pP). He showed (Theorem 3.1(2) in pQ) that 
for F = —X* in R 2n , S(u) has locally finite n-dimensional Hausdorff measure if u 
£ C 2 . We obtained the same result as Lemma 5.4 in |3] by a different argument 
(we used only elementary linear algebra and the implicit function theorem in the 
proof; also we were not aware of pQ at the time was written) . In this paper, we 
generalize this result to the situation of general F (see Theorem D below and its 
proof in Section 6). For u <E C 1 ' 1 and F = —X* in R 2n , Balogh showed (Theorem 
3.1(1) in pQ) that dirriES(u) < 2n — S where diiriE denotes the Hausdorff dimension 
with respect to the Euclidean metric and 5 depends on the Lipschitz constant of 
Vu. He also proved the existence of u S no< Q <iC 1,Q! such that S{u) has positive 
Lebesgue measure for any F £ C 1 (0) where C R m is a given bounded domain 
(Theorem 4.1(2) in pQ). In this paper, we consider functions u of class W 1 ' so that 
the size of S(u) may be large according to Balogh. Therefore for the case of q = 1 
in H1.4|) . we can not neglect the contribution of the singular set to define the weak 
solutions (see Definition 3.2) to the Euler-Lagrange equation of T q : 

(1.5) div VU+ / = H. 

\Vu + F\ 2 -i 

Equation l|1.5|l has been studied in various situations. For F = 0, H = 0, (|1.5|) 
is known to be the q-harmonic equation for 1 < q < oo while it is the equation 
associated to the least gradient problem for q = 1 (see, for instance, ^HJ> |12| . |llj . 
etc.). Geometrically there is a dichotomy for the 1-form O = dz + Fjdxj associated 
to the vector field Vu + F, where F = (Fj). Namely, the hyperplane distribution 
defined by the kernel of might be either integrable or (completely) nonintegrable 
(0 is called a contact form in this case). When F = 0, this is the integrable case. 
For the nonintegrable case (e.g. F = —X*), the quantity of the left side in (|1.5|) 
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with q = 1 can be realized as the p-mean curvature of the graph defined by u in 
pseudohermitian geometry (see l|2.12|l ). 
We study equation (|1.5|) with q = 1: 



(1.6) div Vu + F = H 

\Vu + F\ 

Definition 1.1. Let 51 be a domain in R m , m > 1. We say u £ C 2 (0) is a C 2 
smooth solution to l|1.6[) if and only if l|1.6fl holds in il\S(u). 



In |3j and [2], we considered C 2 -smooth solutions u to (|1.1|) (i.e., (|1.6|) with 
F = —X*) with H = in dimension 2 and, among other things, we proved a 
Bernstein-type theorem. Later in [5] the authors obtained a similar Bernstein- 
type theorem through a different approach. The description of the singular set for 
a C 2 -smooth solution to (|l.ll) occupies a central position in As a geometric 
application, we can show the nonexistence of C 2 -smooth, closed surfaces of genus 
> 2 with bounded p-mean curvature in any pseudohermitian 3-manifold. In [3] we 
also proved a uniqueness theorem for C 2 -smooth solutions for the Dirichlet problem 
of 11.6JI in R 2n . Recently Ritore and Rosales proved a rigidity result for C 2 -smooth 
surfaces of nonzero constant p-mean curvature and an Alexandrov-type theorem in 
the 3-dimensional Heisenberg group (see Theorem 6.1 and Theorem 6.10 in |18| . 
respectively). 

In this paper we consider W ' minimizers for 



(1.7) T{u)= I {|Vu + F| + Hu]dx\ A dx-i A ... A dx m 

f (|1.3l) with q = 1). In Section 3 we define and show that in the space W ' , a 
minimizer for (|1.7|) is a weak solution to the equation i|1.6fl and vice versa (see 
Theorem 3.3). In order to overcome the trouble caused by singular sets which are 
not negligible, we introduce the notion of "regular value". Suppose u € W 1 ' 1 , ip € 
W ' . Define u e = u + eip for e 6 R. We prove that there are at most countably 
many e's for which 

JS(u e ) 

(cf. H1.4(l for q = 1). We call such an e singular, otherwise regular. That is, 
the above integral vanishes for almost all (regular) e (see Lemma 3.1). So we do 
not need to worry about the size of the singular set for regular e's. The idea of 
considering regular values plays a central role both in the proof of the equivalence 
between minimizers and weak solutions and in the proof of the uniqueness theorems 
in Section 5. 

In Section 4 we prove the existence of a Lipschitz continuous minimizer for .F(-) 
with a given boundary value in the case of H = under the following condition on 
F: 



(1.8) 



d K F I =d I f K , I,K = l,...,m 
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for C 1 -smooth functions /if 's (cf. H4.11J) ). We require to be a p-convex domain 
(see Definition 4.1). 

Theorem A. Let Q be a p-convex bounded domain in R m ,m > 2, with dfl S 
C 2 ' Q (0 < a < 1). Let ip S C 2 '"(fi). Suppose F £ C 1 -"^) satisfies the condition 
J _?.£)) for j /or C 1 ' 01 -smooth and bounded Jk's in fl. Then there exists a 

Lipschitz continuous minimizer u £ C 0,1 (0.) for with H = such that u = tp 
on d£l. 

We note that a C 2 -smooth bounded domain with positively curved (positive 
principal curvatures) boundary is p-convex. Also condition (|1.8(l includes the case 
F = —X*. We can actually find out all the solutions to H1.8[) (see H4.13J0 . We notice 
that, for n = 1, Pauls f|16|) proved the existence of a continuous W 1,p minimizer 
for X(-) under the assumption that the graph of the prescribed boundary function 
ip satisfies the bounded slope condition (see [§])■ 

The idea of the proof of Theorem A is to invoke Theorem 1 1 .8 in [§] for a family of 
elliptic approximating equations (see also ^H])- Namely we first solve the Dirichlct 
problem for the following equations: 

V« + F 

Q e u = div( — ) = in Q, 

\Je 2 + \Wu + F\ 2 

u = ip on dfl 

(see Q4.1| ). We end up obtaining a uniform C 1 bound for solutions to the above 
equations, and a subsequence of solutions converges to a Lipschitz continuous min- 
imizer as e — > 0. In Section 4 we give the details of the proof. 

In Section 5 we tackle the problem of uniqueness of minimizers in the Heisenberg 
group of arbitrary dimension (see Theorem B). We also generalize the comparison 
principle in 3 (cf. Theorem C, Theorem C there) to a weak version and for a wide 
class of F's (see Theorem C below). 

Theorem B. Let fl be a bounded domain in R 2n . Let u,v £ W 1,2 (fl) be two 
minimizers for T{ ) such that u — v £ W ' Suppose H £ L°°(f2) and F 

£ W 1,2 (fl) satisfying divF* > (a.e.). Then u = v in fl (a.e.). 

We remark that in the specific case F — —X* the assumptions in Theorem B 
are satisfied. On the other hand, the condition divF* > is essential in Theorem 
B. Let f2 = B2 — B\ C R 2 where B r denotes the open ball of radius r. Consider 
the case F = and H = -. Let u = f(r), v = g(r), and f ^ g with the properties 
that /(l) = g(l), /(2) = g(2), and /' > 0,g' > for 1 < r < 2. Then it is easy to 
see that u and v are two minimizers for the associated (see also page 162 in 
0). 

Theorem C. Let fl be a bounded domain in R 2n . Let F (a vector field) £ 
W 1,2 (fl) satisfy divF* > (a.e.). Suppose u,v £ VF 1,2 (f2) satisfy the following 
conditions: 



MINIMIZERS FOR THE P-AREA 



5 



divN(u) > divN(v) in Q, (in the weak sense); 
u < v on dQ. 

Then u < v in Q. 

In Sections 6, we study the relation between C 2 -smooth solutions and minimizers. 
In (Theorem B there), we proved that if u is a C 2 -smooth solution to (|1 . 1|> in 
dimension 2 with H bounded near a singular point po, then cither pq is isolated 
in S(u) or there exists a small neighborhood B of po which intersects with S(u) 
in exactly a C 1 -smooth curve T through po (the condition on H can be weaker). 
Moreover, T divides B into two disjoint nonsingular domains B + and B~ , and 
N(u)(p+) = lim pe 

b+^ P o N(u)(p) and N(u)(p ) = hmp eB -^ po N(u)(p) exist. Also 
iV(w)(pJ) = —N(u)(pq) (see Proposition 3.5 in J3. ). In Section 6 and the first part 
of Section 7 (see Proposition 6.2, Theorem 6.3, l|7.1|l . and (JZ2J), we will generalize 
such a situation and give a criterion for u to be a minimizer. In particular, suppose 
u is C 2 -smooth. Then Proposition 6.2 or Theorem 6.3 gives a criterion for u to be 
a minimizer in the situation H m -i(S(u)) > while if H m -i(S(u)) — 0, u must be 
a minimizer (see Lemma 6.1). 

In J^j, Pauls constructed two different C 2 (in fact C°°) smooth solutions to 
the p-minimal surface equation with H — or (jl.6|l with H — and i 7, = 

— X*) with the same C°°-smooth boundary value and the same p-area in f2 c R 2 ■ 
These two solutions do not satisfy the criterion in Proposition 6.2 or Theorem 6.3, 
hence none of them is a minimizer. We can also see this fact according to Theorem 
B (uniqueness of minimizers). In Section 7 we construct the actual minimizer for 
Pauls' example (see Example 7.3). 

In dimensions higher than 2, the situation is quite different. The size of the 
singular set can be relatively small under a suitable condition on F = (Fi). For 
x > 0, let [x] denote the largest integer less or equal than x. In Section 6 we obtain 
an estimate for the size of the singular set and a condition on F for a C 2 -smooth 
solution to ljl.6|l to be a minimizer (see Theorems D and E below). Recall that 
dirriE denotes the Hausdorff dimension with respect to the Euclidean metric. 

Theorem D. Let fl be a domain in R m . Suppose u £ C 2 (il) and Fj E C 1 (S1). 
Then for any p S Jl, there exists a neighborhood V of p in such that S(u) D V 
is a submanifold of V satisfying 

(1.9) dim E (S(u) n V) < m - [ rank ^jFj -djFj)(p) + 1 ^ 

Theorem E. Let Q be a bounded domain in R m , m > 2. Suppose u 6 C 2 (f2) 
n C°(0) is a C 2 -smooth solution to with H € C°{n\S{u)) n L°°(ft) and Fj 

€ C 1 (r2). Suppose there holds 

(1.10) [ rank (Mi^Mil+l] > 2 

for all p 6 Q. Then u is a weak solution to J^l.b)) and a minimizer for jJ. 7| ) if in 
addition u E W 1 ' 1 ^). 
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Corollary F. Let Q, be a bounded domain in R 2n . Suppose u E C 2 (fl) n C°(f2) 
is a C 2 -smooth solution to the p-minimal surface equation with H = 0). 

Then in dimension > 4 [n > 2), u is a weak solution to the p-minimal surface 
equation and a minimizer for with H = if in addition u G W ' (CI). 

In Section 8 we study the uniqueness of solutions to elliptic approximating equa- 
tions Q e u = H (see Ij4.1|l ^. e > 0. Since this is an elliptic equation for a given e > 0, 
the uniqueness of solutions follows essentially from the known elliptic theory (see 
e.g. [5]). But for the reader's convenience, we include a proof here as the Appendix. 

We were aware of the paper J2j while this work was being done. After this 
paper was submitted, we were informed of the work |18j . Some problems related to 
this paper were studied in JJ] and ^Hj- We are grateful to Andrea Malchiodi for 
many discussions, in particular, in the study of Example 7.3. We would also like 
to thank the referee for stimulating comments and pointing out many grammatical 
errors. 



2. Hypersurfaces in the Heisenberg group 

In this section we introduce some basic notions for a hypersurface in a pseu- 
dohermitian manifold. By viewing the Heisenberg group or R 2n+1 as a suitable 
pseudohermitian manifold, we give geometric interpretations of (1.1) and (1.4). 

Let (M, J, 8) be a (2n + l)-dimensional pseudohermitian manifold with an in- 
tegrable CR structure J and a global contact form 8 such that the bilinear form 
G = |<i8(-, J-) is positive definite on the contact bundle £ = ker8 (EI)- The 
metric G is usually called the Levi metric. Consider a hypersurface E C M. A 
point p G E is called singular if £ coincides with TE at p. Otherwise, p is called 
nonsingular and V = £ n TE is 2n — 1 dimensional in this case. There is a unique 
(up to sign) unit vector N e £ that is perpendicular to V with respect to the Levi 
metric G. We call N the Legendrian normal or the p- normal ("p" stands for " pseu- 
dohermitian" ) . Suppose that E bounds a domain in M. We define the p-area 
2n-form A by computing the first variation [A will be computed below for the case 
of the Heisenberg group), away from the singular set, of the standard volume in 
the p-normal N : 

(2.1) 6 fN f 8 A (dO) n = c(n) f fA 

Jn is 

where / is a C°°-smooth function on E with compact support away from the singular 
points, and c(n) = 2™n! is a normalization constant. The sign of N is determined 
by requiring that A is positive with respect to the induced orientation on E. So we 
can talk about the p-area of E by integrating A over E (which might not be closed 
from now on). Then we define the p-mean curvature H of S as the first variation 
of the p-area in the direction of N : (the support of / now is also assumed to be 
away from the boundary of E) 

(2.2) S fN [ A = - [ fHA. 

Consider the Heisenberg group viewed as a (flat) pseudohermitian manifold 
(R 2n+1 , 6o, Jo)- Here 8o = dz+ Y^=i( x jd%j' — Xjtdxj) at a point (X, z) = 
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(xi, xy, x n , x n >, z) G R 2n+l and Jo{e.j) = ey , Jo(ey) = —ij where 

d d , d d 

(2 - 3) 6j - d*~ + a;j 'd? 6j ' ~ dx~: Xj d~ z 

j = 1,2, ...n span £ = kerOo- Let S be a graph defined by z = u(X). Note that 
e.j 's and ey 's form an orthonormal basis with respect to the Levi metric Go = 
(5Z" = i dxj A dxj>) (•, Jo-). Observe that an element w = X)j=i( a i^i + byey) G 
Co n TS satisfies d(z - m(X)) (v) = 0. It follows that 

n 

(2.4) ~ x i') a i + + = °- 

Let TV ee -D- 1 X)"=i[( u xj ~ a;j')ej + + ^)ej'] where D = {^2] =1 [{u Xj - 
x j') 2 + ( u x -, + Xj) 2 }) 1 ^ 2 - It is easy to see that N is perpendicular to £ n TE by 

(2.4) , that it is of the unit length w.r.t. Go and hence N is the p-normal (that the 
associated A is positive will be shown below). We can now compute 9o A (d<d ) n 
= c(n) dz A dx\ A dxy A ... A dx n A dx n > and 

(2.5) l n {Q A (dQ ) n } = -ci^D-^il) + (II) + (III)} 
where Ln means taking the interior product with N and (dxi deleted) 



(I) = 'y~][(u Xj — Xj')xj> — (u Xjl + Xj)xj]dx\ A dxy A ... A dx n A dx n > 
i=i 

n 

(II) = — ^(« Xj — Xj')dz A dxi A dxy ...dij A dxy ... A dx n A dx n > 

n 

(III) — y~^(M x<) + Xj)dz A dxi A dxy ...dxj A dxy ... A <ix„ A <ix„'. 
It follows that 

(2.6) 6 fN [ e A(de )" 

Jn 

= f L fN {Q A(de ) n }= f d(i fN {Q A(dQ ) n }) 
Jn Jn 

= JjiN{QoA(de r} 

by the formula L v = t v o d + d o t v and Stokes' theorem. Substituting (2.5) into 

(2.6) and comparing (2.6) with (2.1) gives 

(2.7) A = -D- 1 {(I) + (II) + (III)} 

which simplifies to Ddxi Adxy /\...Adx n Adx n > on S (z = u(X)). Next we compute 

(2.8) S fN A= L fN A = / LfN dA. 

Jt, Jn Jn 
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Here we have used Stokes' theorem and the condition that the support of / is away 
from the singular set and the boundary of E. Noting that D = \Vu — X*\ where 
X* = (xi>, —xi, X2>, —X2, x n >, —x n ), we can easily deduce that d(D^ 1 (I)) = 
and 



d{D + (III)]} = (div — — ^— )dz A dxi A dx v A ... A dx n A dx n >- 

|Vm — X*\ 

It follows from (2.7) and (2.5) that 

(2.9) i N o dA = -(div V " - — )A. 

\Wu-X*\ 

Substituting (2.9) into (2.8) and comparing (2.8) with (2.2), we obtain the following 
expression for the p-mean curvature H-% of the graph £ = {{X , u(X))} : 

(2.10) Hv = div- V " V 



\Vu-X*\ 

Next we consider a general vector field F = (F[) instead of —X*. Let Op = dz 
+ Fjdxi where I ranges over 1, 1', n, n' . Assume that Qp is a contact form, 
i.e., Op A (dO p) n ^ everywhere (satisfied for F = —X* as shown previously). 
For instance, the condition is equivalent to dFy/dxi — dFi/dxy ^ in the case 
n — 1. Define 

d d 

(2.11) ei ■ = — ■- Fi—, I = 1,1',. ..,n,n'. 

It is easy to see that annihilates the e/'s. Define the CR structure Jp on the 
contact bundle kerO^ by Jp{ej) = ey and Jp(ej>) = — ej for j = 1, 2, n. For 
the 2-dimensional case (n — 1), we can find a nonvanishing scalar function A (= 
2(dFy /dxi — dFi/dxy)^ 1 ) such that {ei, ey} forms an orthonormal basis with 
respect to the Levi metric Gp associated to (Jp, XOp). Let ip = z — u(x\,xy) 
be a defining function for the graph of u. By a formula in Section 2 of [3], we can 
compute the p-mean curvature Hp with respect to the pseudohermitian structure 
(Jp, XOp) as follows: 



(2.12) Hp = -div, 



d u xx +F x d Uxy +F V 
dxi \\7u + F\ dx v |Vu + F| 

,. Vu + F 
= div- 



\Wu + F\ 



Here we have used |V b V|G- = \/(eiip) 2 + (eyip) 2 = \Vu + F\ by flTtf) . 
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3. MlNIMIZERS IN THE HEISENBERG GROUP 

In this section we deduce some properties of a minimizer in the Heisenberg group. 
In fact we consider a more general area functional (this is just (11.71) ): 

(3.1) F(u)= [ {\Vu + F\+Hu} 

Jn 

where ft C R m is a bounded domain, F is an arbitrary (say, L 1 ) vector field on fi, 
and H £ L°°(Ti) (we omit the Euclidean volume element). 

Definition 3.1. u £ W 1,1 ^) is called a minimizer for F(it) = J n {\Wu + F| + 
Hu} if T{u) < T{u + ip) for any tp £ Wq' 1 ^), where F G and H G L°°(fi). 

We are going to investigate the first variation of T. Let u,ip € W 1 ' 1 ^) and 
m £ = it + etp for e G F. It follows that u £ — u~ e = (e — tt)ip. Let S(u e ), the singular 
set of u e , denote the set of points where Vu £ + F = 0. So from (|3.1|) (noting that 
| Vu £ + F | = | e - 1 1 I V<p I on 5(u g )) we have 

(3.2) ^(« e ) = \e-e\ f | V<p | + / | + F | 

Js(u f ) Jn\s(u t ) 

+ / i?ttg+ / (e-e)Hip. 
Jn Jn 

Since | Vu e + F | 2 — | Vug + F | 2 = 2(e - e)(Vug + F) • V99 + (e - e) 2 | Vtp | 2 , we 
compute from l|3.2|) 



T{u E ) - FQxg) = [g-el /■ , V | + /" 2(V Mg + F)-Vy+( £ -£)|Vy| 2 
- : " : £ -£ Jsiu,) Ja\s(u t ) \Vu e + F\ + \Vu£ + F\ 

+ [ Hep. 
Jn 

Note that the integrand of the middle term in the right-hand side of the above 
formula actually equals (e — it)~ 1 (| Vu £ + F \ — | Vug + F |) whose absolute value 
is less than or equal to | V<£ | . Therefore by Lebesque's dominated convergence 
theorem, we can easily take the limit as e — > e± (+: the right-hand limit; — : the 
left-hand limit), and obtain 

(3.3) dH j' £±) =±f \V<p\+ I N(u- e )-Vp+ [ Hip 

ae JS{u t ) Jn\s(ui) Jn 

where N(v) = ^+| is defined on Q\S(v). Note that N(u s ) ■ Vtp G L^fiV^ug)) 

since \N(u^) ■ V<p\ < |iV(ug)| |V</?| = |V^| and V<p £ L 1 (fi) by the assumption. 
Also from the above argument, we have the estimate 

■^)-^)l < / \V V \ + \\H\U I 
\s-e\ J n J n 

Namely, T{u e ) is Lipschitz continuous in e for cp £ W ^(Q). Let n(e) denote the 
Lebesque measure of the set S(u e ) n {V(p ^ 0}. We claim that there are at most 
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countably many e's with k(e) > for a fixed ip. First observe that S(u ei ) fl S(u £2 ) 
C {V^ = 0}, and hence (S(u El ) n {V(^ ^ 0}) n (S(u £2 ) n {V^ ^ 0}) = (empty). 
Let |f2| denote the volume of the bounded domain f2. So the number of e such that 
k(s) > i for any positive integer is at most [fi|fi|] + 1 where [x] denotes the largest 
integer less than or equal to x. Therefore there are at most countably many e's with 
k(e) > 0. We call such an e singular, otherwise regular (i.e., k(e) = 0). By l|3.3J) . 
we obtain (|3.4|1 in the following Lemma. 

Lemma 3.1. (1) T{u e ) is Lipschitz continuous in e for ip € W 1 ' 1 ^). (2) There 
are at most countably many singular e's. (3) For a regular e, f g , u ■> \ Vip |= 0, 

d:F ^ u ^ exists, and 

(3.4) <mp± = [ N{Ue) . Vtp+ f Hi p. 

Next for £2,£i regular with £2 > £1, we compute the difference of ^J" 8 ^ for 
e = £2, £1 by (|3.4|) . Using k(sj) = 0, j = 1,2 to shrink the domain of the integral, 
we obtain 

(3.5) ^ - ^1 = / _ jv (u £1 )] • V* > 0. 

Here we have used Lemma 5.1' (also holds for u,v € W 1, ) in 3 to conclude the 
last inequality in [|3.5JI by noting that V<£> = (£2 — £i) _1 (Vii E2 — Vu El ). We have 
the following result. 

Lemma 3.2. (1) rf "^" e ^ is an increasing function of e for e regular. (2) Let Ej, 
j = 1,2,..., be a sequence of decreasing (increasing, respectively) regular numbers 
tending to e (e may be singular) as j — > 00. Then we have 

(3.6) hm — — = ( = , respectively). 

j^oo de de de 

Note that we have the precise expressions for the right-hand limit and 
the left-hand limit ^"^"""^ at e in (|3.3|l . 

Proof. (1) follows from l|3.5|l . To prove (2), first observe that f s , u ■. \V<p\ = 
by the definition of £j being regular. Therefore we have 



(3.7) / |Vy>|=0. 

Ju=^ 1 S(« ej ) 

Let 5*00 = U°^ 1 S*(u £j ). Since \N(u Ej )\ < 1, we estimate | J s ^N(u Ej ) ■ Vtp \< 
f \\7tp\ = by l|3~7jl . So we obtain 



(3.8) / N(u Ej )-Vp = 0. 

It then follows from l|3.4|) and l|3.8jl that 
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(3.9) = I -V(</.,)- V y - - / ll r - 



n\s(u Sj ) 



N(u £ .)-Vtp+ / Hip. 
o\S3o Jn 

On the other hand, observe that lim^oo N(u e .) = iV(ug) in f2\[5oo US(ug)] and 
(3.10) iV(, £ .)= i^ + FJ + ^-e)^ = fi-g)VfP 

in S(u£)\Soo. Now we compute 



(3.11) / iV( W£ .)-V^ 



(/ +/ )N(u Ej )-V<p 



\ £ 3 _ £ l Js(u~ E )\s x JaxiSooUSim)] 



|Vp|+/ N(u ej )-Vip 



-> ± / |V^| + / tf(u 8 ) • Vy> 

JS(ui) Jn\s(ui) 

as j — * oo (+ for decreasing e^; — for increasing e_y). Here we have used <|3.10l) and 
Lebesque's dominated convergence theorem. By (|3.9|) . (|3.11|) . and in view of l|3.3|l . 
we have proved l|3.6|) . 

Q.E.D. 

Definition 3.2. Let C R m be a bounded domain. Let F be an Lj oc vector 
field on f2. Let H G L ; 1 oc (0). We say u € Vl^ 1 (SI) is a weak solution to the equation 
(|1.6|l . i.e., divN(u) = H in 17 if and only if for any ip S Co°(f2), there holds 

(3.12) / |V^|+ / N(u)-V(p+ [ Hcp>0. 

Js(u) Jn\s(u) Jn 



Recall that N(u) = ^+| , fiT(«) = {Vw + i? = 0}, and iV(u) ■ e L x (0\5(u)) 
since |jV(«)-Vy| < |JV(w)| \V<p\ = \V(p\ and V<y3 e L 1 ^) by assumption.. Note that 
with <p replaced by — <p in (|3.12l) . we also have — J s i u \ |Vy| + Jn\s(«) ' 
+ Jo iJ<p < 0. Moreover, if the (m — l)-dimensional Hausdorff measure of S(u) 
vanishes, then the equality holds in (|3.12|l . We remark that in Definition 3.2 for the 
case H € L°°(f2), the space Co°(f2) of test functions can be replaced by W ' (f2) 
since the former is dense in the latter in the W 1 ' 1 norm ( 9 ). Note that in the 
definition of a minimizer, we require u £ W 1 ' 1 ^), F € ^(fl), and H £ L°°{Vt) 
while for the definition of a weak solution, u can be in a larger space W 1 ^), F S 
Lj oc (n), and H G 1^(0). 
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Theorem 3.3. Let u e W 1 - 1 ^), F € L\n), and H € L°°(£l). TTien wisa 
minimizer for F(-) if and only if u is a weak solution to the equation divN(u) 
II . 

Proof. Suppose u is a minimizer for T{u). Then ^£2±i > q, and hence (ETT21 
follows from (|3.3[1 (letting e = in I!OJl). So it is a weak solution. Conversely, 
suppose u is a weak solution. Since T{u e ) is Lipschitz continuous in s by Lemma 3.1 
(1), d:F ^^ exists a.e. (in fact at least for regular e) and it is integrable. Moreover, 
we have 

(3.13) T{u + y) - T{u) = f 1 de. 

Jq de 

On the other hand, from Lemma 3.2 and the definition of weak solution (Defini- 
tion 3.2), we obtain that d "^" e - ) > for any regular e G [0, 1] in view of l|3.3|l (take 
e = 0). By Lemma 3.1 (2), > a.e.. It follows from (|3~T3|) that F(u + ip) > 

T(u). That is to say, u is a minimizer for J-(u). 

Q.E.D. 



4. Existence of minimizers-proof of Theorem A 

Let f2 be a bounded domain in R m ,m > 2. Consider the following elliptic 
approximation u = u e (e > 0) (a geometric interpretation can be found in [lfip 
with given boundary value <p (G C 2 ' Q (f2), < a < 1, say) : 

Vm + F 

(4.1) Q e u = div( ) = in Q, 

'e 2 + \Wu + F\ 2 



u = f> on <9f2 

where F = (Fj), 7=1, ...,to. In the case of m = 2n, I ranges over 1, 1', n, n' 
(e.g., Fj = —xii for the case of a p- minimal surface. Here we use the convention 
that Xj» = —Xj, j = 1, n). We will make use of Theorem 11.8 in [H] to solve 
(14.1(1 in C 2,a (f2) (then a subsequence of u £ will converge to what we want). First we 
check that Q e is elliptic. A direct computation shows that (summation convention 
applies) 



(4.2) 



u H {e 2 + \Vu + F\ 2 ) - [uj + F/)(uj + Fj)ujj 

U E u = - 

[e 2 + |Vu + F| 2 ] 3 / 2 

t (e 2 + [Vu + F\*)diF! - {ui + F!)(uj + F / )9 J F J 

[e 2 + |Vu + F| 2 ] 3 / 2 

= au(s i x, Vu)uu + 6(e, x, Vu) 

where 

, n + [Vu + Fj 2 ) - (u f + F f )(nj + Fj) 

(4.3) a u (e,x, Vu) = 



[£ 2 + |Vm + F| 2 ] 3 / 2 



and 
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„ v (e 2 + \Vu + F\ 2 )d I F I -{u I + F I ){uj + Fj)d I F J 
Ms.x.vu) — - . 

[e 2 + |Vu + F| 2 ]3/2 

For ^ (pi) e R™, we compute from <|4.3I) that 

(e 2 + |Vu + F\ 2 )pj - («j + Fj)(uj + Fj) PiPj 



(4.4) a u p lP j 



[e 2 + |Vu + F| 2 ] 3 / 2 
[e 2 + |Vu + F| 2 ] 3 / 2 



9 9 

>0 



Here we have used Cauchy's inequality |(Vu + F) ■ (pi)\ 2 < |Vu + F\ 2 p 2 (noting 
that p 2 means the sum E/p 2 ). It follows from (|4.2|) and l|4.4|) that Q E is elliptic. 

To apply Theorem 11.8 in we need to get an apriori estimate in C 1 (f2) -norm 
at least. Suppose u e is a C 2 ' a (D.) solution to the equation Q e u = (assuming 
F € C l ' a ((l); later replacing F by a\F), u = aip on < a < 1. In the case of 
Fi = ~x r , djF! = 0, (m + F!){uj + FrfdjFj = -(uj> + Fj,){uj + Fj) = 0, and 
hence b(e, x, Vu) = 0. Since Q £ is elliptic, it follows from the maximum principle 
(see e.g. Problem 10.1 in 0) that 



(4.5) SUp lligr I < Sup |lt e | = SUp \<7<p\ < SUp \<f\. 

n an an an 

Note that the right hand side is independent of e. For a general F, we will invoke 
the comparison principle for a second order, quasilinear operator with a "tail" term 
(namely, Theorem 10.1 in 9 ). First we can find the comparison functions as shown 
below. Let || ||oo denote the supremum norm. Let Br denote the ball of radius R, 
centered at the origin. 

Lemma 4.1. Let C Br C R m be a bounded domain. Suppose F <E C 1 (fi) be 
such that Fj and diFj are all bounded in f2. Then there are C°° -smooth functions 

w _ p x 1 +KR. _|_ e X2+KR w ' _ R. _ e x 2 + K R. 

in R m , where k = «(e, R, HF/Hoo, | \diFj\ \ 00 ) > and n' = k'(e, R, ||F/||oo, 
1 1 9r-Fj- 1 j cso ) > 0, such that Q e w > and Q e w' < in fl. Moreover, we can choose k 
and k' independent of e (but depending on Eq) for < e < So, a positive constant. 

Proof. Let w have the above expression with k to be determined later. Let w\ 
= d Xl w, Wu = d 2 w, W12 = d X2 d Xl w, and so on. It follows that 

(4.6) ton = W! = e Xl+KR , w 22 = w 2 = e X2+KR , and 
wu = 0, otherwise. 

In view of l|4. 2\ with u replaced by w, we compute the dominating (will be clear 
soon) term in the numerator, which is cubic in w as follows: 
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(4.7) w2 j)C^2 w u) - ^2 WiWjWij 

J I I, J 

— w\ w 22 + w 2 Wn (by (|4.6|) ~) 

_ e 2x 1 +x 2 +3KR _|_ e 2x 2 +x 1 +3KR E^J) 

It is easy to see that any other term in the expansion of the numerator is bounded 
by either ae 2KR , c 2 e nR or c 3 for n large. Here Cj = Ci(e, R, \\FiWaa, | |<9/Fj| |oo), 
i = 1, 2, 3, are independent of k. Therefore we have Q e w > in by 1)4. 7|l for a large 
K = k(e, R, \\FiWoo, \\diFjWoo). Moreover, n is independent of e for < e < £q, 
a positive constant. Similarly, we can find n' = k'{e, R, ||i*/||oo, ||9/-Fj||oo) >0 
such that Q e w' < in fi. 

Q.E.D. 

Proposition 4.2. Let ft C Br C R m be a bounded domain. Let F G C 1 (fi) 
such that Fx and diFj are all bounded in fi. Suppose u £ G C 2 (fl) PI C (f2) satisfies 
14- If) , i.e., Q e u e — in 0, and u e = o~(p G C°, < a < 1, on <9f2. TTien i/iere exists 
a constant C = C(e, i?, ||-Fj||oo, ||<9/-Fj||oo, IMIoo) (independent of a) such that 

(4.8) sup|u £ |<C. 

n 

Moreover, the bounds hold uniformly for < e < Sq, a positive constant. 

Proof. Let w, w' be the comparison functions as in Lemma 4.1. On dVL, w < 
o-Lp + Ci = u e + d for some constant C\ = Ci(e, R, H-F/Hoo, H^jFjUoo, IMIoo) 
(independent of a) independent of e for < e < £q> a positive constant. On the 
other hand, we have Q e w > = Q £ {u e + C\) in ft by Lemma 4.1 and the observation 
that Q e (u e + C\) = Q e u £ . Now we apply the comparison principle for quasilinear 
operators (e.g. Theorem 10.1 in 9 ) to conclude that 

(4.9) w < u £ + d in fl. 

Similarly, there is a constant C 2 = C 2 (e, R, ||-F/||oo, | |<9/i<j| |oo, IMIoo) (inde- 
pendent of a) independent of e for < e < £o, such that w' > a<p — C 2 = u £ — C 2 on 
dfl and Q e ui' < = Q e (ue — C 2 ) in CI. So we obtain from the comparison principle 
that 

(4.10) w' > u e - C 2 in n. 
Thus (O follows from g^J) and (OU- 

Q.E.D. 

For the gradient estimate, we will reduce the problem to a gradient estimate at 
the boundary. We need to require a condition on F. Suppose there are C 1 -smooth 
functions /r-'s (K = 1, m) in CI such that 



(4.11) 



d K Fi = dif K . 
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We remark that if both Fi and Gi satisfy the condition l|4.11(l . so does Fx + Gj. 
In fact, we can write down all the (local) solutions to (|4.11(l . It is easy to see from 
gJH that d K (djFi - diFj) = for all I, J,K = 1, m. It follows that 



(4.12) 



djFx - diFj = Cu 



where the constants Cxj satisfy the skew- symmetric relation: C/j = —Cjx- Since 
the left-hand side of 1(4.12(1 is linear in F, the general solutions are the solutions to 
djFj — diFj — plus a special solution. Let ui = J^i Fjdxi. Then du> — if djFj 
— diFj = 0. So locally there is a function g such that u> = dg. Hence Fi = dig. On 
the other hand, we observe that Fx = \ YIk Cik x k is a special solution to l|4.12|l . 
So the general solutions to (|4.12|l are 



(4.13) 



Fi = dig + Ci K x K - 



K 



It is then easy to verify that F = (Fi) having the form 14.13JI are also solutions to 
JUID for f K = d K g + \ TjCjkXj- 



Proposition 4.3. Let ft C R m be a bounded domain. Let F — (Fi) £ C 1 (fi) 
satisfy the condition \4-.ll\j , for all I,K — 1, m, where all fxc's are bounded. 
Suppose u e G C 2 (f2) satisfies the equation Q s u £ = Hq, a constant, in fl. Then we 
have 



(4.14) 



sup \d K u e \ < sup \d K u e \ + 2\\f K \ 
n an 



Proof. Write Vu + F = (m + Fi). Let D £ (u) = \J e 2 + |V« + F\ 2 . Compute 
(summing over J while fixing / and K) 



(4.15) 



d 



ui + Fi 



K- 



D s {u) 

Uik + d K Fi (ui + Fi)(uj + Fj)(u JK + d K Fj) 



D e {u) 
5 U - vi{u)vj{u) 
D e (u) 



dj{u K + fx) 



where vi(u) = (ui + Fi)/D 6 (u) and we have used the condition 14.11(1 . Now for 
v £ Cq(O), we compute 
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(4.16) = / (Q £ u £ - H )d K v = [ di ^ Us J ] I , + ^* d K v (summing over /) 

Jn Jn D £ {u e ) 

= - / did K v 
Jn D e (u £ ) 

= [ d K^rr~T^ d i v ( d i d K = d Kdi) 

Jn De{u £ ) 

{ajj(e 7 x,'Vu £ )dj[(u £ )K + fic]}9l v (summing over / and J) 

by Ij4.15|l with u replaced by u £ . Here au(e, x, Vm e ) = [Su — v i(u £ )v j(u £ )]/ ' D £ {u £ ) 
(cf. Ij4.3|) ). It is then easy to see that (|4.16(l holds also for v € Co(^) (use the 
regularization Vh of (7.13) in [3j to approximate v). So (u £ )k + /k is a weak 
solution to the equation Lw = di{ajj(e,x,\7u £ )djw} = (cf. (8.2) in [5]). By 
(I4.4|i . this is an elliptic equation in divergence form. So by the maximum principle 
(e.g. Theorem 8.1 in |5] with b l = c 1 = d = and ajj bounded), we have 

sup \{u £ ) K + Ik\< sup \{u £ ) K + fal- 
si an 

Then (|4"T4) follows. 

Q.E.D. 

For a general F, the bound for Vu e may depend on e if we invoke the maximum 
principle for a more general situation (for instance, Theorem 8.16 in 9 ). 

To perform the boundary gradient estimate, we need a comparison function to 
apply the comparison principle. Let fl C R m be a bounded domain with coordinates 
denoted by x\, X2, x m . We call a coordinate system orthonormal if it is obtained 
by a translation and a rotation from x\, X2, x m . We define a certain notion of 
convexity for £1 as follows. 

Definition 4.1. We call f2 C R m parabolically convex or p-convex in short if for 
any p £ dfl, there exists an orthonormal coordinate system (ii, £2, i m ) with 
the origin at p and fl C {ax\ — X2 < 0} where a > is independent of p. 

Note that a C 2 -smooth bounded domain with positively curved (positive princi- 
pal curvatures) boundary is p-convex. 

Proposition 4.4. Let 51 C R m be a p-convex bounded domain. Suppose u £ <E 
C 2 (f2) n C 1 (H) satisfies Q £ u £ = in Q and u £ = ap G C 2 (Cl) on d£l with F £ 
C 1 (f2) for < a < 1. Then there exists a constant C = C(e, a, ||-Fj||oo) H^/FjUoq, 
||<9/</?||oo, I \didjip\ |oo) (independent of a) such that 

(4.17) sup|Vu e |<C*. 

an 

Moreover, the bounds hold uniformly for < e < £0, a positive constant. 

Proof. Given p € <9f2, we have an orthonormal coordinate system (x%, X2, x m ) 
as in the definition of p-convexity. Consider the comparison function w = aG + otp 
where G is the function G = ax\ — X2 viewed as a function of (2/), 1=1, 2, to, 
for large a to be determined. In view of the invariance of Q £ (u) under the coordinate 
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changes of translations and rotations, we compute (Q E , D £ being the corresponding 
operator, quantity of Q £l D e with respect to (xj), respectively) 

(4.18) Qe( w ) — Qe(u>) (w is w viewed as a function of (xj)) 

P(G)a 3 + Aa 2 + Ba + E 
D £ (wf 

by (|4.2|) where P{G) is the corresponding quantity of P(G) = G Xi G X2X2 — 2G Xl G X2 G XlX2 
+ G X2 G XlXl with respect to (xj), and A is a function of a, Fi, diFj, dip, didjip 
while B, E are functions of e, a, Fj, diFj, diip, didjip. Moreover, a direct com- 
putation shows that P(G) = 2a. Since a > 0, Q E {w) > (<, respectively) = Q E (u e ) 
for positive (neg ative , respectively) large a = a(e, a, ||-Fj||oo, ||9/-Fj||oo, ||0j¥>||oo, 
||9/9j^||oo) by (Qj) . Note that a is independent of cr and independent of e for 
< e < £q. On the other hand, w = aG + atp < (>, respectively) atp — u £ on 9ft 
since G < on ft by the p-convexity. Therefore w < (>, respectively)it e in ft by 
the comparison principle for second order quasilinear operators (see e.g. Theorem 
10.1 in 9 ). Noting that G(p) = and hence w(p) — crip(p) — u E (p), we then have 

4.19 = < (> respectively — 

ov ov 

where v = —d X2 at p. Observe that ^ (for either positive or negative a) is bounded 
by a constant depending on s, a, H-P/Hoo, H^i-FjIIoo, H^Hoo, \\didj(p\\oo, but 
independent of cr, p (moreover, the bounds hold for < e < £q, & positive constant), 
so is by 14.19JI . Since u e = op on 9ft and u e — op S C 1 (0), we can easily 
show that the derivatives of u £ — op in the xi, X3, x m ) (except $2) directions 
all vanish at p. It follows that in the Xj (j 7^ 2) direction, the derivative of u £ is 
the same as the derivative of op. So of course it is bounded by ||Vy>||oo (note that 
< a < 1). Altogether we have proved H4.17JI . 

Q.E.D. 

Proof of Theorem A. 

In order to apply Theorem 11.8 in [Hj to solve the Dirichlet problem 1)4. l[l . we 
consider a family of equations: 

(4.20) Q £j(T u ee div ; V " + = in ft 

+ |Vii + (rF| 2 



u = atp on 9ft, < o < 1. 

Express Q e . a u — a/j(e, 1, V«; cr)uu + 6(e, X, Vu; cr) where au(e, x, Vit; a) and 
o(e, x, Vu; cr) are given by l|4.3|> with F replaced by aF. It is then easy to check 
that the conditions (i), (ii), (iii) on page 287 of |2] are satisfied. To have an apriori 
Holder estimate for Vit, we invoke Theorem 13.2 in 9;. Comparing l|4.20[) with 
(13.2) in gives 

A(x, u, Vu) = — ; ^ , B(x, u, Vu) = 0. 

e 2 + \Vu + oF\ 2 
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Following pages 319-320 of [H], we find a IJ = D Pj A ! = au(e, x, Vu; a) and A(e, x, u, 

Vu) = e 2 /[e 2 + |Vw + crF| 2 ] 3 / 2 by Therefore we can take X K = e 2 /[e 2 + (K + 

C)2]3/2 in ^ 13 4 j of p in which K = n (- see page 53 in jgj for the notat i on ) 

and C = ||-F||oo- Similarly we estimate 

\DpjA 1 ] = |a w (e,x,Vu;£r)| < 1 = < -. 

^/e 2 + |Vu + crF| 2 6 

So we can take Ak = e _1 . Since both D Z A T and £? vanish, we compute 



\8 J A I \ + \B\ = \D XJ A'\ 

\{e 2 + \Vu + ( rF| 2 )9 J (a^) - ( Ul + + ^jgjVgQj 

[e 2 + |Vu + (rF| 2 ] 3 / 2 

^e 2 + |Vu + crF| 2 

Therefore we can take an upper bound [i K = e _1 (| + n)supjj |9j-Fj|. Now by 
Theorem 13.2 in we have an apriori Holder bound for Vit in terms of n, K (= 
Mi ; fi), Ak/Xk, HkI^k, size of J7, and |v|2;0- On the other hand, we observe that 
Lemma 4.1, Propositions 4.2-4.4 still hold for Q e , a instead of Q e . So we have an 
apriori C 1 bound for solutions of (|4.2Q(I . independent of a and e (for < e < 6q). 
Altogether we have obtained an apriori C 1 " 3 (ri) (f3 > 0) bound for solutions of 
(|4.2U|I . independent of a (but depend on e). By Theorem 11.8 in [SJ, we obtain 

Theorem 4.5. Let f2 be a p-convex bounded domain in R m ,m > 2, with dQ G 
C 2 ' Q (0 < a < 1). Let ip € C 2 a {Ci). Suppose F € C 1,a (fi) satisfies the condition 
14-11]) for C ,a -smooth and bounded Jk's in fl. Then there exists a solution u £ £ 
C ' a (n) of the Dirichlet problem: Q e (u) = in Cl, u = cp on dfl for given s > 0. 

(Proof of Theorem A Continued) 

Propositions 4.2-4.4 tell us that there exists a constant C — C(e, a, R, H-FjHoo, 
||3/^||oo, \\didjipWoo, ||//||oo) such that 



(4.21) sup |u e | + sup |Vu e | < C. 

Moreover, the bounds hold uniformly for < e < £o, a positive constant. 

In view of (|4.21() we can find a subsequence u e . (0 < < eq, £j — > 0) converging 
to uo in C° by the Arzela-Ascoli theorem. Then the Lipschitzianity of uq follows 
by taking the limit of ratios: (x =/= y) 



u Ej (x)-u Ej (y) 
x-y 

Next we claim that uq is a minimizer for (see l jl.3^ 1 such that uq = (p on 
dft. Observe that W 1,q (Vt) is compactly imbedded in L 1 (fi) (e.g., Theorem 7.26 in 
[Hj). So we may as well assume that u £j converges to uq in L 1 (J7). Also note that 
\p + F\ is convex in p since \Xpi + (1 - X)p 2 + F\ = |A(pi + F) + (1 - A)(p* 2 + ^)| 
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< X\pi+F\ + (1 - X)\p 2 +F\ for < A < 1. We can therefore apply Theorem 4.1.2 
in to conclude the lower semicontinuity of (see {OJ0 : 

(4.22) F(u ) <lim inf T{u s \ 
Now for v £ W 1 ' 1 with v — (p G Wq 1 , we estimate 

(4.23) J-(u Ej ) = / | Vu ej + F | (omitting volume element) 

J Jn 

< J y/e* + \Vu Bi +P\* 



< / a/£, 2 + |Vi; + F| 2 
Jo 



< £j uoZ(fi) + / \Vv + F\ 
Jn 

where we have used the fact that the Dirichlet solution u e - £ C 2 (Vl) is also a 



mimmizer for T €] (u) = J n ^e] + |Vu + F\ 2 . Taking the limit infimum of l|4~2*3*|l 

and making use of H4.22|l . we finally obtain that F(uq) < J n \Vv + F\ = T(v). That 
is to say, uq is a minimizer for F(-). 

Q.E.D. 



5. Uniqueness of minimizers-proof of Theorems B and C 

Recall (see Section 3) that ft C R m denotes a bounded domain and T(u) = 
$ n {\Vu + F\ + Hu} for u £ W 1 ' 1 ^), F G L 1 ^), and H G L°°(fi). We will prove 
two (W 1 ' 1 ) minimizers for .F(u) with the same "boundary value" have the same 
normal vector " almostly" . 

Theorem 5.1. Let u, v G W ' (fi) 6e two minimizers for T(u) such that u — v 
£ W Q ' (CI). Let u e = u + e(v — u). Then for any pair of regular e±, £2 G [0,1], there 
holds N(u El ) = N(u 62 ) in Q\[S(u El )U S(u E2 )](a.e.). 

Proof. By (|3.13|) with <p = v — u, we have 

(5.1) = T{v) - T(u) = [ ^P±a!e. 

Jo de 

As in the proof of Theorem 3.3, the same argument shows that '^"^ > for 
any regular e £ [0, 1]. In view of jnij and Lemma 3.2(1), = for any regular 

e £ [0, 1]. It follows from that / n \ s( „ e) N(u e ) ■ V(v - u) = 0. Therefore for 
any pair of regular ex, £2 G [0, 1], there holds 



(5.2) f [N(u E2 )-N(u El )}-V(v~u) = 0. 

Jn\[s(u ei )us(u e2 )] 
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Here we have used f„, s N c , , N(u £o ) ■ V(u — u) = and f„, c , > N(u £ , ) ■ 

J S(u E1 )\S(u e2 ) \ £2J \ ) J S(u C2 )\S(u cl ) V £ i' 

V(« - u) = by observing that for j = 1,2, |iV(u Sj ) • V(w - u)| < \V(v - u)\ 
and f S ( u \ |V(u — n)| = from the definition of £j being regular. Write v — u = 

(u £2 — u ei )/(e2 — £i) for e 2 ^ e±. By Lemma 5.1' in [3], the integrand in (|5.2|l is 

|Vu £2 + F| + |Vu ei + FL . , , ... , |2 



■\N{u e2 )-N{u ei )Y 



2(e 3 -Ei) 

It then follows that N(u £l ) = N(u £2 ) in f2\[£(tt El ) U «S(ttg 2 )]. 

Q.E.D. 

For a vector field G = (.91,32, •••,<?2n) on C i? 2 ™, we recall that G* = (32, —31, 

34, -33, 32n, -02n-l). 

Lemma 5.2. Let u, i> 6 where the domain ft is contained in R 2n . Let 

u £ = u + e(v — u). Suppose N(u ei ) = N(u £2 ) in Q\[S(u £l ) U S(u £2 )} for a pair e%, 
£2 such that £1 ^ £2. Then for j — 1,2, there holds 

(5.3) (Vu £] + F)* ■ (Vv - Vm) = in (a.e.). 

Proof. We will prove l|5.3[) only for j = 1 (similar argument works also for j = 2) 
For p e S(u £l ), Vu £l + F = 0. So (|S"3|) holds obviously. For p £ S{u £2 ), (j5~3|) also 
holds by observing that Vw-Vu= [(Vu El +F)- (Vu E2 + F)]/(ei -e 2 ) = (Vu £l +F 
)/(ei — £2) and G* ■ G = 0. For the remaining case: p £ fi\[S(u £l ) U 5(u e2 )], we 
observe that for j = 1,2, 

(5.4) N(u £] y ■ Vu Ej = |- = F* • N(u ej ). 

\Vu £j + F| 

Here we have used the property G* • G = twice. Since N(u £l ) — N(u £2 ) in 
f2\[5(u El ) U S(u £2 )] by assumption (hence N(u £l )* = N(u £2 )* also), we take the 
difference of (15 .4|) for j — 1 and j = 2 to obtain 

(5.5) W(u ei )* • (Vu E2 - V« ei ) = 0. 

Formula l|5.3|) for j = 1 on £l\[5(w £l ) U S(u £2 )] then follows from 1)5. 5|) by noting 
that v — u = (u £2 — u £l )/(£2 — ex). 

Q.E.D. 

We will use the following general criterion to prove the uniqueness of minimizers 
and a comparison principle for weak functions later. 

Theorem 5.3. Let Q, be a bounded domain in R 2n . Let w £ W ' p (Cl), a £ 
W ,9 {fl), where 1 < p < 00, q = (q = 00 for p — 1). Let F (a vector 
field) £ W 1 ' 1 (Q) D L g (Q.) satisfying divF* > (a.e.) or divF* < (a.e.). Suppose 
(Vcr + F)* ■ Vw = in (a.e.). Then w = in Q (a.e.). 
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Proof. Take Uj € C§°(Q) -> w in W 1 ^ and F k e C°°(fi) ^ F in VF 1 ' 1 n L q . 
Suppose ujj does not vanish identically. Then there exists a decreasing sequence 
of positive numbers a, converging to such that = {\u)j\ > a.;} CC fl is not 
empty for large i and dSljj is C°° smooth (by Sard's theorem; note that \ojj\ is C°° 
smooth where uij ^ 0). Also we take v k €E C°°(Vl) — > cr in W 1,2 . Consider 



(5-6) hi,k,k= / W (V^ + Ffc)*-^ 

where v denotes the boundary normal. We first compute 



(5.7) / |^-| {Vv k + F- k )* -v = a, (Vv k + F k )*-v 

JdQ jti JdQ 



div[(Vv k )* + F£] 



= at I divF^ 

Here we have used Green's theorem for the second equality and div(\7v k )* = 
for the third equality in lj5.7Jl . It follows from (|5.7(l that 

(5.8) lim I. a s = 0. 

I— »oo 

On the other hand, a similar reasoning gives 



(5.9) I M = / VK-|-(V^ + F s r + |^| d»v[(V« fc r+J^] 



Observe that 11;%,.* = > 0} = Q\{oJj = 0}, (ft\{ w j = 0})\a, Vi = Ug i (%,j+i 
S W 1 ' 1 ^), and hence 



(5.10) ( / ){V| Wi | • (V« fe + + 

if!,,, in\{uj=o} 

= / {V|w 3 -| • (V« fc + F s )* + | WJ | di«#|} 

by ijEBJ. It follows from l(53ty. (f^TTPfc that 



o = / v|^-| • (y Vk + P- k y + luji divPi 

JQ,\{ujj=Q} 

I VK-|-(V^ + ^r + |^-| divF-*. 

JQ. 



Here we have used V|wy| = if u>j =0 (p. 152 in |2]). Letting k — > oo in the above 
formula gives 
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(5.11) = / V|w 3 -|-(V«k + F)* + |w i | divF*. 

Letting k — > oo in the first term of l|5.11|l . we then estimate by using the as- 
sumption Vw ■ (Vcr + F)* = 

(5.12) f V\lUj\ ■ (Va + F)* 
J si 

/ (Vu>j - Vw) • (V(T + F)* - / (Vu>j - Vw) • (Vcr + F)* 

J{u)j>0} J{u)j<0} 

— > as j — > oo. 

Here we have used LOj — > w in and (Vcr + F)* S L 9 (r2) by assumption. For 

the second term of (|5.11|l . we have 

(5.13) lim / \wj\ divF* = / |w| divF* > or < 

if w (noting that divF* > or < by assumption). By l|5.11|l . I|5.12|l . and 
(|5.13|l . we reach a contradiction. Therefore w = in f2 (a.e.). 

Q.E.D. 



Remark. If F does not satisfy the condition in Theorem 5.3, then the theorem 
may not hold as shown by the following examples. Let £1 = (0,7r) x (0,7r) C R 2 . 
Let w = sin x sin y € W . Then Vw = (cos x sin y, sin a; cosy). Take a — and 
F = (cos a; sin y, sin a; cosy). It is easy to see that F* = (sin a; cosy, — cos x sin y), 
divF* = 0, and (Vcr + F)* ■ Vw = F* ■ Vw = 0. With the same <r (= 0) and w as 
above, we can also take F = sin x (cos x sin y, sin a; cos y). Then still (Vcr + F)* ■ Vw 
= F* ■ Vw = while divF* = cos a; sin a; cosy has no definite sign in f2. 

Proof of Theorem B. 

The proof follows from Theorem 5.1, Lemma 5.2, and Theorem 5.3 with p = q = 



2, cr = u Sl , and w = v — u. 



Q.E.D. 



Next we want to prove a comparison principle for weak sub- and super- solu- 
tions (a comparison principle for C 2 -smooth functions has been studied in . See 
Theorem C and Theorem C there). First we need to define relevant differential 
inequalities in some weak sense. Let il C R m denote a bounded domain. Recall 
that N(u) = ^ v "+^ is defined on Cl\S(u) (F, say, is an L\ oc vector field in f2). 

Definition 5.1. Let H £ L\ oc (Vl). We say u S W 1 ' 1 ^) satisfies divN(u) > H 
(< H, respectively) in the weak sense in f2 if and only if for any ip G C^°(fl) and 
ip > 0, there holds 
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(5.14) - f \Vip\ + [ N(u)-V<p + 
Js(u) Jn\s(u) 



Hp < 

(5.15) (/ \Vcp\ + I N(u)-Wp + I H(f > 0, respectively). 

Js(u) Jn\s(u) Jn 

Recall that we defined the weak solution to divN(u) = H in Section 3 (see 
(|3.12|l ). The following result justifies the above definitions. 

Proposition 5.4. Let H e L°°(Cl). Then u S W X ' X {Q) satisfies divN{u) > H 
and divN(u) < H in the weak sense if and only if u € W 1,l (Cl) is a weak solution 
to the equation divN(u) = H. 

Proof. Since Cg°(fi) is dense in Wq' X (SI), (|5~TH) and (|5~T5j) hold for every 
(p £ C^°(Cl) if and only if they hold for every cp G Wq 1 ' 1 ^). Write ip = p + — ip~ for 
ip £ W ' (Cl) where p + = max{^>, 0} and p~ = max{— ip, 0}. Express 



(5.16) f \Vip\+ [ N(u)-\7p + [ Hip 
Js(u) Jn\s(u) Jn 

= {[ |V^+|+ / N(u)-Vp+ + [ Hp+} 
Js(u) Jn\s(u) Jn 

-{- f \V<p~\+ [ N(u)-Vp-+ f Hp-}. 
Js{u) Jn\s(u) Jn 

Note that ip + > and p~ > 0. Now suppose u is a weak solution to divN(u) 
< H and divN(u) > H. Then the right-hand side of (|5.16l) is nonnegative by 
our definitions. So the left hand side of (|5.16l) is nonnegative, i.e., (|3.12|) holds. 
Conversely, suppose u is a weak solution to divN{u) — H. That is to say, the left 
hand side of (|5.16|l is nonnegative (note that p is not restricted to be nonnegative 
here). By taking ip > i.e. p~ = (<p < i.e. p + — 0, respectively) in l|5.16|) . we 
obtain (|5.15|) f (|5.14|) . respectively). 

Q.E.D. 

Definition 5.2. u,v £ W /1 (J7) satisfy divN(u) > divN{v) in Cl in the weak 
sense if and only if for any p £ W ' (Cl) and p > 0, there holds 

(5.17) -/ \Vp\+ [ N(u)-Vp<+ I \Vp\+ I N(v)-Vp. 

Js(u) Jn\s(u) Js(v) Jn\s(v) 

Definition 5.3. u, v S W ' (f2) satisfy u < v on dCl if and only if (u — v) + = 
max(w- v,0) e W^ifl). 

Theorem 5.5. Suppose u,v £ W ' (Cl) satisfy the following conditions: 

divN(u) > divN(v) in CI (in the weak sense); 
u < v on dCl. 
Then N{u) = N(v) on {u > v}\[S(u) U S(v)]. 
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Proof. Let ip = (u — v) + . The condition u < v on dfl implies that ip <E W ' (ft). 

dj r (v £ 
' de 

111 T omina *-i 

de — de 



Let v £ = v + ep. From Lemma 3.2 (1), rf ^J E ^ is increasing in regular e. It follows 
that djr ( t,Q +) < j^gkJ L emma 3 2 (2). In view of the formula l|3.3[) . we have 



(5.18) + [ \Vp\+ [ N(v)-Vp<- [ \Vtp\+ [ N(m)-V(p. 

Js(v) Jn\s(v) Js(vt) Jn\s(vi) 

Observe that Vi — u on {u > v} and tp = on {u < v}. So the right hand side 
of l|5.18|l equals the left hand side of H5.17JI . It follows that 

(5.19) -/ \V<p\+ [ N(u)-V<p = + [ \V<p\+ [ N(v)-V<p. 

Js(u) Jn\s(u) Js(v) Jq\s(v) 



Write 



(5.20) / N(u)-V(p- [ N(v)-Vp 

Jn\s(u) Jn\s(v) 

(N(u)-N(v)) ■ Vip 

n\[S(u)us{v)] 



N(u) ■ Vip - / N(v) ■ Vip. 

S(v)\S(u) Js{u)\S(v) 



We claim 



(5.21) -/ \V<p\- [ N(v)-Vp = 0. 

Js(u) Js(u)\S(v) 

Since p = on {u < v}, we only have to discuss the case that u > v. In this 
case, ip = u — v and hence Vip = (Vu + F) — (Vv + F) = — (Vv + F) in S(u) 
(and = in S(u) n S(v)). So N(v) ■ Vp = ^+| • [-(Vv + F)} = -\Vv + F\ in 
S (u)\S (v) . It is now clear that (|5.21(l holds. Similarly there also holds 



(5.22) -/ \V(p\+ N{u)-Vtp = 0. 

JS(v) Js(v)\S(u) 

Combining JOty . (|S~2T)|) . jSH gives 



(5.23) / (iV(u) - N(v)) ■ Vip = 0. 
Jn\[s(u)us(v)] 

By Lemma 5.1' in 3 (which works also for u, v <G W rl ' 1 (0)), we have 

(5.24) (N(u) - N(v)) ■ Vp = ^ U + F ^^ V + F ^ N(u) - N(v)\ 2 

on {u > v}\[S(u) U S(v)] (where ip = u — v). Noting that ip = on {u < v} and 
substituting into l|5.23|) . we finally obtain N(u) = N(v) on {u > v}\[S(u) U 

S(u)]. 

Q.E.D. 
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We can now prove the comparison principle for weak sub- and super- solutions. 
Proof of Theorem C. 

By Theorem 5.5 and Lemma 5.2 (switching the roles of u and v and taking il = 
{u > v}, £i = 0, e 2 = 1), we obtain (Vv + F)* ■ V(u - v) + = 0. Then we apply 
Theorem 5.3 (with p = q = 2, a = v, and w = (u — v) + ) to conclude that (u — v) + 
= in. f2. That is to say, u < v in Q. 

Q.E.D. 



6. When a smooth solution is a minimizer 

In this section we determine when a smooth solution is a minimizer. We will 
prove Theorem D, Theorem E, and Corollary F. We first prove a result for the case 
H m -i(S(u)) = 0, in which a C 2 -smooth solution must be a weak solution. 

Lemma 6.1. Let Q, be a bounded domain in R m . Suppose u £ C 1 (r2) n 
C 2 (Q\S(u)) n C°(Cl) satisfies JHJ) in Q\S{u) with F e C 1 (Cl\S(u)) and H e 
C°(fl\S(u)) n L} (Q). Suppose H m ^i(S(u)), the m — 1 dimensional Hausdorff 
measure of S(u), vanishes. Then u is a weak solution to i)j.6|) and a minimizer for 
fTfl if u e W 1 ' 1 ^) and H G L°°{VL) also. 



Proof. By Theorem 3.3, it suffices to prove that for any ip S C^°(il) (3.12) 
holds. That is, 

/ |V^|+ / N(u)-\7(p+ j H<p>0. 
Js(u) Jn\s(u) Jn 

Write SI = fl + U Q U 0_ where S7+ = {(f > 0}, il_ = {ip < 0}, and fl = {ip 
= 0}. If £1+ 7^ 0, then there exists a sequence of ej > approaching 0, such that 
Sl ej = {(p > Ej} 7^ 0, U°^ 1 il £j = il+ and 9f2 £j . are C°°-smooth by Sard's theorem. 
Since u € C 1 (f2), 5(11) n 0, £j is compact. Together with the condition H m _i(S(u)) 
— 0, for any a > 0, we can find a finite cover of balls B rk (pi~) of center p^ and 
radius r^, k = 1, 2, for S*(u) n f2 £j . such that 

(6.1) 5^JZ m _ 1 (SB f . |k (p fc )) <a. 

k=l 

On the other hand we compute by the divergence theorem and the equation (|1.6|) 



(6.2) / (<p- ej )N(u)-v 

Jd(n ej \uB rk ( Pk )) 

Vcp ■ N(u) + (<p - ej)H. 

n Sj \uB rk ( Pk ) 

Since tp — Sj — on dfl e , we can estimate the boundary term in 1)6.2(1 as follows: 
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(6.3) | / (<p- ej )N(u)- V \ 

Jd(n S] \uB rk ( Pk )) 

< {max^ - e j \}H m -i(uf^ 1 dB rk (p k )) 

a max \tp — eA 

Q J 

by $jin\ and the fact that \N(u)\ = \v\ = 1. Letting a — » in gives 



(6.4) / V<p-N(u) + (tp-£ 3 )H = 

Jn e .\s(u) 

in view of Letting ej -> in we obtain 



(6.5) / V<^-iV(w)+ / tpH = 

Jn + \s(u) Jn + 

by noting that the volume of {0 < ip < Sj} tends to as £j — > 0. Similarly we also 
have 



(6.6) / \7ip-N(u)+ [ 



tpH = 0. 

n_\s?(u) Jfi- 

On the other hand, it is obvious that the integral of ipH over Slo vanishes since ip 
= on fig. Observing that V<p = a.e. on Oo in view of Lemma 7.7 in we 
conclude that 

(6.7) / X7tp ■ N(u) = 0. 
It now follows from JfTSJl, JESJI, and (|6~7|l that 

(6.8) / V^-7V(u) + / <pH = 

for ip £ Cg°(f2). Comparing (|6.8|l with (|3.12|) and noting that the first integral of 
(|3.12|l is zero by H m _i(S(u)) = 0, we have completed the proof. 

Q.E.D. 



Proof of Theorem D. 

Write Vu + F = (ui + Fj)^. Consider the map G : p e Q ^ ((«J + -F/)(p))/Li- 
Computing the differential dG of G at a singular point p (where G(p) = 0), we obtain 
(djui + djFi) in matrix form (note that G € C 1 ). From elementary linear algebra 
we compute 



(6.9) 



rank (djui + djFi) + rank (diuj + diFj) 
> rank {{djuj + djFj) - (dm.; + OjFj)} 
= rank {djFj — diFj). 
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Observing that rank (djui + djFi) = rank (diuj + diFj) (the transpose has the 
same rank), we can deduce from (|6.9[) that rank dG(p) > [ rank ^ 2 "(p» +1 ] where 
hji = (djFj - diFj). It follows that 

(6.10) dl m(Ker dG{p)) <m- ["*"* + \ 

Then by the implicit function theorem there exists an open neighborhood V of p 
in £1 such that G _1 (0) n V = S(u) n V is a submanifold of V, having (Euclidean) 
dimension diniE bounded by the right side of (|6.10|1 . 

Q.E.D. 

Proof of Theorem E. 

It suffices to prove that H m _i(S(u)) — in view of Lemma 6.1. Combining (|1.9|) 
and (|X. XOp . we bound dimES(u) by m — 2. It follows that H m _ 1 (S(u)) = 0. 

Q.E.D. 

Proof of Corollary F. 

For m = 2n, F — —X*, we compute rank {hji) = In. Therefore ((1. 10(1 is 
reduced to n > 2, hence m > 4. 

Q.E.D. 

We remark that the condition (|1.10|) does not hold in dimension m = 2. So 
Hi(S(u)) may not vanish. Therefore a C 2 -smooth solution may not be a mini- 
mizer in this case (see Example 7.4). We will discuss the general situation that 
H m -i(S(u)) > below. 

First we will give a criterion for, in particular, a C 2 -smooth solution to be a 
minimizer. Let Q be a domain in R m . Let rcflbeam- 1 dimensional, orientable, 
C 1 -smooth submanifold. Let B CC ft be an open neighborhood of a point in T 
with C 1 -smooth boundary and B being compact. Suppose TtlB divides B into two 
disjoint parts (note that T may or may not contain some singular points). That is, 
B\T = B\ (rnB) = B + UB~~ where B + and B~ are disjoint domains (proper open 
and connected) (see Figure 1(a) or Figure 1(b) below). Suppose u is C 2 -smooth in 
f2\r and has no singular points in fi\r. Let F G C 1 (r2) for simplicity. Suppose also 
N + (u) and N~(u) (restrictions of N(u) to B + and B~ , respectively) are continuous 
up to rn B, i.e., N+(u) e C°(B+), N~(u) € C°(B~-), so that divN*^) = H in 
B ± , respectively. Let v + and v~ denote the outward unit normals to T n B with 
respect to B + and B~ , respectively. Note that = — 

Proposition 6.2. Suppose we have the situation described above. Then u is 
a weak solution to \1.6}) on B with H € C (B\T) fl L°°(B) if and only if along 
rnB, there holds 

(6.11) (N + M - N~(u)) ■ v + = (N + (u) - N~(u)) ■ v~= 0. 

Note that for u S lF 1:1 (i?), u is a weak solution to (|1.6|) if and only if u is a 
minimizer for (|1.7fl in view of Theorem 3.3. 

Proof. Using the divergence theorem, we compute 
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(6.12) 




/ ip(N + (u) - N~ (u j) ■ v 
JrnB 



Here we have used v = —v + and divN(u) = H in both B + and B . Observing 
that H m (S(u) n B) = since H m (S(u) n B) < H m (T (IB) = 0, we conclude from 
(|3.12|l (also ip replaced by —ip) that u is a weak solution to (|1 . f>f> if and only if 

(6.13) / N(u) ■ V(f + H(p = 

Jb\t 

for all ip £ C^(B). On the other hand, (f^TSfl holds if and only if j5TT|) holds by 
(|6~T^ . 

Q.E.D. 

In order to have a criterion for a more general situation, we extend Proposition 
6.2 as follows. Let ft C i? m be a bounded domain. Let A C P C such that T is 
relatively closed in 17, H m _i(A) — 0, and r\A is a C 1 -smooth to — 1 dimensional 
manifold. Suppose f2\r = U^^fl,-, the union of at most countably many domains 
fij. For each j, we have 90j C d£l U T. We can view f2\T as domains flj obtained 
by cutting apart along T and T\A as the union of two copies of T\A. Let Vj denote 
the outward unit normal to d£lj. Then i/j exists for any point p £ dflj D (F\^4). 
At p, there is another I (I may equal j) such that v\ = — Vj. Let F £ C 1 (fi\F) 
for simplicity and iJ € C°(il\F) n L ; 1 oc (Sl). Suppose u £ C 1 (fi\r) has no singular 
points in il\T. Let Nj(u) denote the restriction of N(u) on flj. 

Theorem 6.3. Suppose we have the situation described above. Furthermore, 
suppose Nj(u) £ C°(flj U (dtij n (T\A))) n C 1 ^) satisfies divNj(u) = H in Qj 
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for any j. Then u is a weak solution to J^l.b}) in ft if and only if for each p € r\^4 7 
there exist j, I as described above, such that at p, there holds 

(N^u) - #,(«)) • v.= (N.(u) - #,(«)) ■ v x = 0. 

We should remind the reader that for u £ and H S L°°(f2), u is a weak 

solution to l|1.6|) if and only if u is a minimizer for (|1.7|l in view of Theorem 3.3. 

Proof. Let U CC CI have compact closure in i7, and suppose that the boundary 
dU is C 1 -smooth. For ip 6 Co°(n) with support contained in t/, we compute 



(6.14) / N(u) ■ Vip + Hip = / N(u)-V<p + Hip 

'U JU\T 

(pN(u) ■ v 

d{u\r) 

= I <P(Nj(u) ■ Vj + JV,(u) • u{) 

{ ^ l) Jan J n(r\A)nu 

- E / ipiNjM-Ntiu)).^. 
{ ^ l) Jan J n(r\A)nu 

For the last equality we have used vi — —Vj. Now observe that u is a weak solution 
in £1 if and only if the first term of Q6.14|l vanishes for any p> £ (CI) and associated 
U. On the other hand, this is equivalent to concluding that (Nj(u) — Ni(u)) ■ Vj — 
on T\A bv 

Q.E.D. 

We remark that it is possible that u £ C 1 \C 2 while N(u) £ C 1 in the nonsingular 
domain. For instance, let u = xy + g(y) with g £ C 1 \C 2 . Take F = —X*. We can 
then compute N(u) = (0, ±1) in the nonsingular domain defined by 2x + g'(y) =/= 0. 

We will also make a remark on deducing the second equality in 1|6.14)1 . First 
note that at points of A with H m -i(A) — 0, v may not exist. How do we deal 
with this? For any e > 0, we can find a finite open cover U^ =1 Dj D A such that 

Y^ k j=\ H m -i(dDj) < e. By the divergence theorem we have 



N(u) -Vcp + Hcp= / ipN(u) ■ v. 

l(U\r)\u* =1 Dj Jd[(u\r)\u^ =1 D 3 ] 

Passing to the limit as e — > and observing that the integrands are bounded (since 
|JV(te) | = 1), we obtain 



N(u) -Vip + Hip= / <pN(u) ■ v. 
>u\r Jd(u\r) 

The idea of the above argument was used in 0] . We have displayed this idea in the 

proof of Lemma 6.1. We also used a similar argument in the proof of Theorem 5.2 

in PJ. We remark that Pauls had a similar result (for m = 2, F = —X*, and H 

= 0) as Theorem C in 17]. Ritore and Rosales also obtained a similar result for 

C 2 -smooth minimizers (for m = 2, F = —X*, and H = constant) as Theorem 4.15 

in Q2]. 
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7. Examples 

We shall give examples of Lipschitz (continuous) minimizers in dimension 2. 

Definition 7.1. A p-area minimizer or a p- minimizer in short is a minimizer for 
(Ol) with H = 0. 

Throughout this section, we will always work on the situation that m = 2, 
F = —X*, and H = 0. Recall that the integral curves of A r± (u) are straight lines 
(see Section 4 in [2]), called the characteristic lines, segments, or rays. We call 
the angle between T (oriented) and a characteristic ray (with direction N- 1 ^)) in 
B + (B~ , respectively) touching a point p G T the incident (reflected, respectively) 
angle at p. Therefore geometrically <|6.1 f li is equivalent to saying that at p S TDB, 
either N + (u) = N~(u) (see Figure 1(b)) or N + (u) ^ N~(u) which implies 

(7.1) The incident angle=The reflected angle. 

(see Figure 1(a)). Suppose u G C 2 at a point peTflJ] and L is a singular curve. 
Recall that if the characteristic line segments T + and T_ in B + and B~ respectively 
meet at p, then r + U {p} U T_ must form a straight line segment according to (the 
proof of) Proposition 3.5 in 0. Therefore by (|7.1() (note that N + (u) = —N~(u) 
at p in this situation) , we can conclude that 

(7.2) r + and T_ are perpendicular to T at p if u G C 2 at p. 

The constraint (|7.2|l gives a necessary and sufficient condition for a C 2 -smooth 
solution of Ijl.lD with H = to be a p-minimizcr. We can have a function u G 
C 2 (f2) which satisfies the p-minimal surface equation divN{u) = in fl\S(u), but 
is not a weak solution or a p-minimizer. 

Example 7.1. Consider F — (—y,x) in the following N(u)'s. 

(a) By taking a — cost?, b — sind, and g(~bx + ay) = (cot-d) (—bx + ay) 2 in 

(1.2) of [3] for < ■& < ^, we obtain u(x,y) — — xy + y 2 cot This is a C 2 smooth 

solution to divN(u) = in R 2 \S(u) for F = (— y 7 x) by a direct computation. We 
can easily determine the singular set S(u) = {u x — y = 0, u y + x = 0} = {y = 0}. 
On the other hand, N (u) = (cos - !?, sin??) which is not perpendicular to the a;-axis 
{y = 0} (see Figure 2(a)). So in view of l|7.2f> . this u is not a p-minimizer on any 
bounded domain f2 containing part of the s-axis. 

(b) Let u(x, y) — —xy + y 2 cot # for y > 0; = — xy + y 2 cot r\ for y < 0; = for 
y = where < i9, 77 < 27r, ■& ^ n, 77 7^ n. We compute 

(7.3) Ar x (u) = (^L-|smtf|, Isintfl) for y > 0; 

smw 

N-'-Ou) = ( — C ° S ^ I sin 77 1 , — I sin 77 1 ) for y < 0. 
sin 77 

Observe that (|7.1I) (or 1)6. H[> ^ holds if and only if i? + 77 = 2ir (see Figure 2(b)) by 
(17.3(1 . Therefore we conclude that u is a (C 1 * 1 -smooth) p-minimizer on any bounded 
domain in R 2 if and only if t9 + 77 = 2-7T in view of l|7.1|l . 

Example 7.2. Let u(x,y) = xy for y > 0, and u = for y < 0. Consider the 
case of F = (— y, x). Compute 
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Figure 2. 



N^iu) = (1,0) for x > 0,y > 0; A^u) = (-1,0) for x < 0,y > 0. 
^(«) = - 7 %ll=for y <0 

(see Figure 3). Observe that the positive y-axis {x = 0, y > 0} is a singular 
curve where l|7.2|l holds true. Also on the x-axis {y = 0} except the origin, A rJ -(u) 
is continuous and hence (|6.11(l holds true (note that the x-axis is not a singular 
curve, but is a curve where u is not C 1 smooth). Applying Theorem 6.3 with T = 
{x = 0, y > 0} U {y = 0}, we conclude that u is a (Lipschitz) p-minimizer on any 
bounded domain £1 C R 2 ■ 




Figure 3. 



We remark that it is not possible to construct a Lipschitz p-minimizer having 
a loop consisting of characteristic lines (see Figure 4 for an example). Indeed, by 
contradiction, suppose that the loop consists of three characteristic lines 7 l5 7 2 , 
and 7 3 as indicated in Figure 4. Let A denote the region surrounded by 7 l5 7 2 , and 
7 3 . We integrate the contact form 8 = du + xdy — ydx over the loop as follows: 
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/A 

This contradiction confirms our claim 



/ 9 (7i,72j and 7 3 being Legendrian) 

"'7iU7 2 U7 3 

/ dO (Stokes' Theorem) 
J A 

2 dx A dy = 2 Area(A) ^ 0. 
Ja 



Figure 4. 



Example 7.3. There can be two distinct C 2 -smooth p-minimal graphs (i.e., 
satisfying on nonsingular domain) having the same boundary value and the 
same p-area, but both of them are not p-minimizers. Consider u — x 2 + xy, v = 
xy + 1 — y 2 (first given in ^B]). We can easily verify that u and v satisfy Hl.lfl with 
H = on their respective nonsingular domains and have the same value on the 
unit circle in the xy - plane. But they do not satisfy (|7.2(1 . So by Proposition 6.2 
or Theorem 6.3, neither of them can be a p-minimizer. Compute the p-area (see 
H1.2B 1 of u and v over the unit disc A as follows: 

X{u) = f VS\x\dxdy = 

J a 3 

8V2 



X(v) = / 2\x- y\dxdy 

J A 



I A 

So they have the same p-area. By the uniqueness of p-minimizers (see Theorem B), 
we also conclude that neither u nor v can be the p-minimizer. 

We are going to describe what the (unique) p-minimizer looks like on A with 
the boundary value (or curve) p(6) = cos 2 6 + cos 9 s'm8 (9 is the standard angle 
parameter for dA). Let (a, (3, 7) be a point of a line segment L c A x R meeting 
the boundary curve with the projection L C A passing through the origin. Suppose 
9' is the angle between the positive £-axis and part of L, lying in the upper half 
plane. Then we have 

(7.4) a = t cos 9', (3 = tsm9' 

7 = cos 2 9' + cos 9' sin 9' 
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for — 1 < t < 1 (note that p(ir + 9) = p(9)). Suppose the contact plane passing 
through (a, (3, 7) intersects the boundary curve p at (cos#, sin#, p(9)). Then we 
have the following relation: 



(7.5) p(9) - p(9') + t sin(0 - 6') = 

by observing that z — 7 + x(y — [3) — y(x — a) = is the equation for such a contact 
plane in R 3 with coordinates (x,y,z). By elementary trigonometry for the above 
specific p, we can reduce (|7.5|) to 



(7.6) ^[sin(20+|)-sin(20' + |)]+Mn(0-0')=O. 

The idea is to choose 9 such that sin (26' + f ) = 0. Then we solve £29 for 
8 (perhaps we have multiple solutions). Keeping L or L associated to 6 as the 
singular set in mind, we connect (a, f3, 7) G L to a point of the boundary curve, 
associated to by a line segment. Since these line segments are Legendrian, their 
union forms a Legendrian ruled surface, hence a p-minimal surface ( 3 ). Moreover, 
if two characteristic lines (i.e., above Legendrian lines projected to the xy-plane) 
meet at a point of L, condition (|7.1() holds. So in this way we can construct the 
p-minimizer by Proposition 6.2 or Theorem 6.3. We give more details below. 

First solving sin(26*' + ?) = gives 9 = (n — j)^ where n is an integer. There 
are two such 9 f, s modulo an integral multiple of tt, namely 9' = |7r and 9 1 = |7r. 
We take 6' = |7T (it turns out that 9' = |7r won't give rise to a p-minimal graph 
in the following argument). So l|7.6|l is reduced to 

(7.7) ^ S[n2 (9-K) = tsm(9-^). 

(note that for 9' = ^ir we have —t instead of t in 17.7fl ). By the double angle 
formula, we deduce from l|7.7|l that 

(7.8) (a) co S (8 - ^tt) = -L sin (0 _ * w ) = . 

The solutions to (6) of 1)7. 8fl are |7r +nir for any integer n, which we ignore. We 
have two solutions 6\, 9^ (modulo an integral multiple of 27r) to (a) of (|7.8|l for a 
given t with the relation 



(7.9) ei -^n=^7r-9 2 . 

When t runs from —1 to 1, 9\ runs from |7r to |7r clockwise while 82 runs from 
— |7r to |7r counterclockwise (See Figure 5 below). 

Denote the line segments between (a, (3) € L (9 1 = |-7r) and the boundary 
point (cos9j, sin#j), j — 1, 2, by Tj, Tf, respectively. Tl and are the xy-plane 
projections of two Legendrian lines , Y"l connecting (a, (3, 7) £ L to (cos 6j, sin 8j, 
P(@j))i J — 1' 2, respectively. We will define a graph u over A, whose restriction to 
the region Q = LU (Uj = i : 2--i<t<iYi) is i U {^j=i^--i<t<i^ t )■ We can parametrize 
Tl , say, in the following form (see (4.9) in[Hj): 
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Figure 5. 



x = s(sin 77(f)) + a(t) 

y = -s(cos??(f))+/3(f) 

z = s [/3(f) sin 77(f) +a(t) cos r)(t)] +-f(t). 

Here 77(f) = § + 2 {t) - 8{t) in which cos(0 2 (f) - |tt) = ^ (see (a) of JEHJ) 

and tan (5(f) = f-^/l — f 2 /2/(l — t 2 /^/2) by elementary plane geometry (we leave 
the details to the reader). On the other hand, (|7.1|l holds along L due to l|7.9|l . 
Therefore u S C 1,1 is a weak solution to (1.1) with H = over the region fL The 
remaining domain A\f2 consists of four small fan-shaped regions (see Figure 5). For 
each of such regions, we can connect two points on the boundary curve, indicated 
by 6 and 9 which are related by 1)7. 6JI with f — 1. Thus we obtain a family of 
Legendrian line segments whose lengths are getting smaller when both 0' and 6 
tend to some critical value (e.g., for the fan-shaped region between and |7r, the 
critical value is jir by solving ^ = 0). These Legendrian line segments form a 
portion of the graph u over A\0. So u is a C 2 -smooth p-minimal graph over A\fi. 
Altogether u S C ' (A) is the (unique) p-minimizer by Theorem 3.3. 



8. Appendix: uniqueness of solutions to (|4.1|) 

The existence of solutions to (|4.1|l is asserted in Theorem 4.5. In this section 
we are going to prove the uniqueness. In fact we can obtain more general results. 
First we define 



(8.1) N e (u) 



Vu + F 



Let a = (e, Vu+F) € Rx R m = R m+1 . Denote \a\ as a. Similarly let /3 = (e, Vv+F) 
and/3= 0\. 
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Lemma 8.1. Let u, v £ W /1 (fi) where fl C R m (m > 1) is an arbitrary domain. 
Then 

(8.2) (N e (u) N £ (v)) ■ (V« - Vv) >^| N e (u) N e (v) | 2 . 

Moreover, the equality holds for e — 0. When e > 0, (N e (u) — N £ (v)) ■ (Vu — Vu) 
= if and only if Vu = Vd. 

Proof. We compute 

(8.3) (N s (u) -N s (v)) ■ (Vu- Vv) 
,Vu + F Vv + F 



) ■ {(Vu + F) - (Vv + F)} 



= { M^i-fe^Q } . {(£ ,v„ + #)-( £ ,v„ + #)} 

= {--5}-{3-/3} = (a + j9)(l-cose) 
a p 

where a • [3 = a/3 cos 9. On the other hand, we can estimate 

(8.4) | N e (u) - N e (v) | 2 

< |--|| 2 =2(l-cos0). 

Now (|8.2|l follows from (|8.3|l and l|8.4|l . Observing that the equality in l|8.4|l holds 
for e = 0, we obtain the equality in (|8.2|) for e — 0. Suppose (N e (u) — N E (vj) ■ (Vu 
— Vu) = 0. By (|8.2(l we have N e (u) — N £ (v). Taking the modulus of this equality 
gives | Vu + F\ = \Vv + F\ if e > 0. It follows that Vu = Vu. 

Q.E.D. 

We remark that the equality in l|8.2|l for e = has been obtained as Lemma 5.1' 
in 0. Recall Q E u = divN E (u) (see (IQl . (jHHJ). Note that for F = 0, £ = 1, Q e u is 
the Riemannian mean curvature of the graph defined by u. In this case, the above 
inequality has been obtained in |14| . and independently. 

Definition 8.1. Let ft C R m be a bounded domain and e > 0. Suppose u, v E 
W 1 ^) and F is measurable. We say Q e u — Q e u > (< 0, respectively) weakly if 
for any <p £ Co(f2), (p > 0, there holds 



(8.5) / (N e (u) ~ N e (v)) ■ Vf < ( > 0, respectively). 

Jn 

Note that N £ (u) and N e (v) are integrable since they are bounded by 1. We have 
the following comparison principle for Q e . 

Theorem 8.2. Let O C R m be a bounded domain and e > 0. Suppose u, v £ 
C 1 (f2) n C°(f2) satisfy Q £ u — Q £ v > (< 0, respectively) weakly and u — v < 
f> 0, respectively) on dQ. Then u ~ v < (>0, respectively) in ft. 



3(5 



JIH-HSIN CHENG, JENN-FANG HWANG, AND PAUL YANG 



Proof. Given a > 0, we choose a function f a G C' 1 (i?) with the property that 
f a = in (—oo, a], f a > 0, and f' a > in (a, oo). Observe that f a {u — v) 6 Cq(17) 
(i.e., / a (u — v) G C 1 and has compact support in 17) by the assumption w — v < 
on <917. It follows from (|8.5|l that 

(8.6) > / (N £ (u) - N E (v)) ■ V(/.(« - v)) 

= f (N s (u)-N s (v))-f' a (u-v)(Vu-Vv) 

J {u—v>a} 

> (by fOjl l. 

Therefore we have (V« — Vu) • (-/V e (u) — N £ (v)) = in {u — v > a} since f' a {u — v) 
> and (N £ (u) - N £ (v)) • (Vit- Vu) > in It follows that Vu = Vw in {u - 

w > a} by Lemma 8.1. Thus we obtain u ~ v = a in {u — v > a}, a contradiction. 
So {it — v > a} is empty. Since a > is arbitrary, we conclude that {u — v > 0} 
is empty. So u — v < in 17. 

Q.E.D. 

We remark that basically the above result can be deduced from Theorem 10.7 
in!. 

Corollary 8.3. Let 17 C R m be a bounded domain. Let e > 0. Suppose u, v G 
C 2 (17) n C°(fi) and F G C^O) safe/?/ Q £ u = Q £ v in 9, and u = v on <917. Then 
u = v in 17. 



Theorem 8.4. Let 17 C i? m 6e a bounded domain and e > 0. Suppose u, v G 
W ' (f2) satisfy Q £ u — Q £ ti > f< 0, respectively) weakly and (u — v) + ((u — 
v)~ , respectively) G W ' (17). XTien it — « < f> 0, respectively) in 17. 

Proof. First we observe that (|8.5jl still holds for tp G Wq 1:1 (S7), tp > 0. It follows 
that 

(8.7) > I (N e (u) -N s {v)) ■ V(it - v) + 

= / (JV e («)-iV e (t;))-V(« - «) 

J{«-u>0} 

> 

by (|8.2|) . Therefore if u — u > 0, V(w — w) + = V(w — f) = by Lemma 7.6 in 

(JH3, and Lemma 8.1. Also if u - v < 0, V(u - u)+ = by Lemma 7.6 in 
Altogether we have shown that V(it — v) + = in 17. Now applying the Sobolev 
inequality to (u — v) + G W ' (17), we obtain (it — t>) + = in 17. That is, u — v 
< in 17. 

Q.E.D. 

We remark that the proof of Theorem 8.4 is based on the idea of the proof of 
Theorem 8.1 in 0. 
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